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1. mscoocncii 


One of the difficult problems in experimental fluid dynemics remains the 
deteiminaticxi of the vort.icity field in fluid flows. The difficulty arises 
from the fact that vorticity is a quantity defined in terms of locad velocity 
gradients 2 md the available velocity meeistirement techniques e.g. Laser Doppler 
Velocimetry and Hot Wire Anerocmetry, 2 ure only sensitive to the local flow 
velocity. Current methods for vorticity measurements use multiple probes to 
ind^jendoitly measure one or more coc^nents of the velocity vector at two 
closely spaced spaticLL locations, from vdiich a vorticity Vcilue is inferred by 
applying some kind of finite difference scheme (1,2). These methods suffer 
fron sane dravbacks. Firstly, vhen Hot-Wire probes are used, the local flow 
vorticity may be strongly influenced due to non-negligible blockage effects 
caused by the multiple probes. Secondly, the locaticxi of each Individual 
mecisuring probe volume (LDV and H(f alike) may not be spatially close enotgh to 
resolve and accurately measure the loccil velocity gradiants. As a consequence 
the measured vorticity is only a "spatially filtered" estimate of the actual 
vorticity field. Another inportant aspect that needs to be considered is that 
these techniques can only provide "one point" information. In order to obtain 
Whole field data, meausuranents must be cairried out sequentially one point at a 
time. Although this sequentiail method can be easily inplemented in 
applications involving steady flows, it is of rather difficult application in 
unsteady flow phenanena. 

Recently, a novel velocity measurement technique, cairoonly known ais Laser 
Speckle or Particle Image Displacement Velocimetry became available (3-7). 

This technique permits the simultaneous visualization of the two-diroensicxial 
streamline pattern in unsteady flcMs and the quantification of the velocity 
field. 

The nain advantage of this new technique is that the whole 
two-dimensional velocity field can be recorded with great accuracy and spatial 
resolution, fron which the instantaneous vorticity field can be easily 
obtained. This constitutes a great asset for the study of a variety of flows 
that evolve stochasticadly in both ^jace and time, and in the case of 
interest, to the study of unsteady vorticeil flews which occur in rotorcraft 
and in high angle of attack aerodynamics. 
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2. KtQKHFIf GF TEE TBCOliaUB 


Hie 2 ^plicaticxi of Laser SJjeckle or Particle Image Displacemeit 
Veloclmetry to fluid flows followed the use of Laser Speckle methods in solid 
mechanics, in vhich c±>ject displacement and rotaticMi were accurately measured 
( 8 ). 

Vtien applied to the measurement of fluid flows the technique involves two 
ec[U2dly inportant steps. In the first step a selected plane of the flow field 
is recorded. Similarly to Laser Doppler Veloclmetry the flow is seeded with 
amad.1 tracer particles. A sheet of coherent light is used to provide a 
"surface" within the fluid. A pulsed Ictser such as a Ruby or a NdYag laser, 
or a CM laser with a Gutter device is nozmally used as the light source. Hie 
laser sheet is formed, for exarple, by focusing the leiser beam first with a 
long focal length spherical l«is, to obtain minimum thickness, euid then 
diverging the beam in one dimaision with a cylindriced. lens (Fig. 1). Hie 
light scattered by the seeding piarticles in the illuninated plane provides a 
mcving pattern. Mhen the seeding concentration is low, the p)attem consists 
of resolved diffraction limited images of the piarticles. When their 
concentration increases, the images overlap and interfere to produce a random 
spieckle pattern. A multiple exposure photograph records this moving pattern. 

In a second step the local fluid velocity is derived from the ratio of 
the meaisured spacing between the images of the seme tracer, or speckle grain, 
eind the time between exposures. Hie recorded image formed by isolated Airy 
disks, in the cases of low seeding ooncaitraticxi, or speckle grains for high 
seeding ooncentration, is a ccnplicated randan pattern. In this image the 
local displacements can hardly be measured by visual or ccttputer-aided 
inspection. Methods for the direct analysis of these images have been used 
but with limited success (9). Hie displacement meeisurements have an 
associated error vhich becomes ixrportant vhen the spacing between successive 
tracer images becomes snail. Hiis method also fails whai the mean distance 
between tracers is of the same order of magnitude as the distance a tracer 
particle travels during the time between exposures. As a consequence the 
trex:er ocxicentration is usually kepjt small, resulting in velocity measurennents 
with poor spatial density. 

Alternative methods used in similar data ctnedysis in solid mechanics are 
applicable and provide good results. Hra methods are currently being 
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developed and investigated, niese are the local coherent illixnination to 
produce Young's fringes and local tvo-dimensional correlation of the image. 

In the first method, the diffraction produced by local oc^erent ill unination 
of the multiple imstges in the photograph generates stredght fringes in the 
Fourier plane of a lens. Ihese fringes have an orientation vhich is 
perpendicular to the direction of the local di^lacenent vector and a spacing 
vhich is inversely proportional to its magnitude (Fig. 2) . Hie displacement 
vector is determined by the aned.ysis of these fringe patterns. Amcxigst the 
various anailysis methods being developed, Fourier and Autooorrelation 
techniques seem to be the most effective and accurate. In the tvo-dimensional 
correlation method a anall area of the image being interrogated is digitized, 
in a NXN point format (vhere N is the nunber of pixel rows or columns) , and a 
two>dimen8ional correlation is performed. This results in a digital 
autocorrelation function with a maxima at the coordinates corresponding to ^ 
average displaceroait of the tracers (Fig. 3). Basically this method is 
equivalent to ^e general Young's fringe method because the two-dimensional 
fourier transform of the fringe pattern is equal to the autooorrelation 
function. The major dravhack of this method is that the computation of the 
autooorrelation function requires large data arrays end beocines extremely slow 
vhen N (nunber of pixel rows or colunns) is large. A new processing method, 
proposed by Yao cuid Mrian (10), reduces the genered NXN element of a 
tvo-dimensional problem to two N elero^it one-dimensional problems, by 
ccn^ressing the information in two orthogoned directions using integratico 
techniques. Although attractive, this procedure leads to seme loss of 
information. Therefore, vhat has to be balaioed is this loss versus the gain 
in computational speed. Future plans in this research include a critical 
evaluation of each of the above mentioned methods. However, we believe that 
the image analysis using the Young's fringe method is the most premising 
approach. 

2.1 Image An 2 dysis Using the Young's Fringe Method 

Consider the function D(x,y) describing the light intensity produced by 
the scatterers in the image pleine of the photographic camera, vhere x,y are 
the plane <x»rdinates . Considering that there is a in-plane displaceroa:it, dy, 
of the scatterers, the speckle or particle images will be translated by M3y, 
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v^ere M is the magnification of the cenera lais. The resulting intensity 
distribution in the image plane beccmes 

D(x,y) + D(x,y4Mdy) = D(x,y) Ot C6(x,y) + 6(x,y-tMdy)] (1) 

vhere 6(x,y) is the Dirac delta function centered on x,y, amd considering that 
a translation can be represented as a ccxwolution vath a delta function. The 
total intensity is recorded on a photographic plate. After development, the 
trananittahce, r , of the negative is given by 

r = art) D(x,y) ® C 6 (x,y + Miy)] (2) 

vhere a euid b are chauracteristic ccxistants of the photographic enulsion. 

TMO-dimensional Fourier transformation of the image is efficiently 
performed by optical processing of the photographic negative. This is 
achieved by locally analyzing the film negative with a prcbe laser beam. This 
process produces, as explained in the following > a straight fringe p)attem 
vhose orioitation and magnitude are related to the directim auid magnitude of 
the displacement. The resulting Fourier transform of the trananittance 
distribution, F , is 

7(u,v) = a 6 (u,v) + b D(u,v) [1 + exp( i2 n vMiy/X )] (3) 

cl 


vhere F represents the Fourier transform of F , u and v are angular 
coordinates of a point in the Fourier plane of the lens, X is the 

wavelength of the interrogating laser light beam. 

The first term, a 6(u,v) on the r.h.s. of equation 3 represents the image 
of a point source, i.e. the interrogative beam, vhen diffraction effects are 
neglected. This image is seen as a smedl bright spot in the center of the 
Fourier plane (Fig. 2). The second term is ccnposed of a diffraction pattern, 
characteristic of the speckle or particle images, modulated by 
Cl+exp(i2wMdy/ X ^ ] 

The intensity distribution for the second term is ctotained by multiplication 
with its conplex conjugate, resulting in 

|0(u,v)|^ C 4 cos^( TTvMdy/x^)] (4) 

a 
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Ihe diffuse background, lC(u,v) I is modulated by a set of Young's fringes 
vAK>se simcing is given by 

X ff 

Where f^ is the fbczd length of the oorveiging lens. Knowing M, X^ and 
measuring dp, the dia^e^zement is easily found fron equation 5, with the 
direction of the motion perpendicular to the oriaitation of the fringes. 

Ey scanning the multiple ejqposed photograph, one can resolve the tvo 
oo np o n ents of the dlsp>lacanait and derive the velocity vector at every point 
in the field. Ihis is a unique capability of this technique. 

2.1 Nature of the Present Researdi 


Ihe motivation for the present vork is to determine and develop the 
necessary background for the use of laser Speckle or Particle Snage 
Disp>lacanent Velocimetry to study vortical unsteady flows . 

In this phase of the work we report on the optimization of the 
photographic parameters involved in obtaining the multiple exposure 
photograph. 


3. A CTTTTKT.nt FOR THB CBDICB (F TEE miOGBAmC EABMdERS llilVXVED IN 
A lASDBBBir 05TH5 lASBt SBBCOZ (LSV^) OR EMTIGI2 SflOB DISEUCaSHT 

VBLucnaiBr (riif) 

The peuTzmeters playing a role in velocity measurements using laser 
Speckle or Particle Snage Displacement Velociinetry are: 

(i) Tracer: type, ccncoitration and dimensions. 

(ii) Exposure parameters: time between e^qxsures, ^posure time and nuiiser 

of exposures. 

(iii) Film parameters: sensitivity, grain and resolution. 

These parameters are interrelated and their choice depends v^xxi several 
factors ais follows: the flow field of interest, namely the velocity range 
and area of the flow to be recorded; the required spatial resolution and 
accuracy; avail^le hardware such as type of laser, CW or pulsed, laiser power 
and laser licpit wavelengidi, and on the camera used for the photography. 

3.1 Tracer 

Tracers to be used in a measuranent involving ISV or PIDV have to be aus 
anall as possible to acxnirately follow the fluid motions, and good light 
scatterers. These requirements aure usually met by tracers used in UN 
applications. In liquid flows we reocmnend the use of pliolite particles, 
which aure practicadly neutradly buoyant in water, latex and titauiiun dioxide. 
In air flows, oil amoka or fog produces a relatively uniform seeding. 

As previously mentioned, the light scattered by the seeding particles can 
create, depending on the seeding concentration, tvo aitirely different 
patterns. This originates two modes of operation of the technique. For lew 
seeding oonoentrations, the pattern consists of resolved diffraction limited 
imaiges of the particles. When this concaitration increases, the images 
overlap and interfere in the image plane to produce a random speckle pattern. 
In either case, the pattern is recorded in a nultipLe exposed photograph fran 
which the velocity data is obtadned. However, in the Paoticle image mode of 
operation, due to light seeding, regions of the flow field may be left 
unseeded, or with poor seeding, resulting in signal drop>-out. Cn the other 
hand, in the Laser Speckle mode of operatixxi, the drop>-out problem is 
minimi but other more restrictive ocnplications may occur, am described 
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belw. Tliese problans precauie the use of the Laser Speckle mode of operation 
in the majority of £Lov/ fields of interest. 

The Laser Speckle mode of c^>eration relies on the recording of identical 
laterally shifted speckle patterns. Slight out-of-plane motion of the 
scatterers, due to flow three-diroensionedity, between the exposures results in 
speckle patterns that are not entirely similar. As a consequence, their 
correlation decreases and local coherent illunination does not produce Young' s 
fringes, or prodvxies fringes with poor contrast and signal to noise ratio 
(SNR) . This poses a severe limitation in the use of the Laser ^jeckle 
Oisplacemoits for the measurement of velocity in turbulent flows, or flows 
with an important velocity ocnqpcxient in the direction perpendicular to the 
measuring plane. The fringe quality becomes less d^sendent on the out-of- 
plane motion vhen individual particle images are imaged and recorded, i.e. in 
the Particle linage mode of operation. In this case the tolerance to 
out-of-plane motioxi is roughly equivedent to the width .of the illumination 
sheet and focusing depth of field of the recording optics. In addition, the 
high concentrati<xis required by the Speckle mode of operation cause multiple 
scattering effects vhich result in the spreeuiing of the thin laser 
illuninating sheet. Furthermore, the high concentration of tracers may 
strongly influence the flow field. 

In conclusion, for a successful measurement using this technique, there 
are {practical bounds to seeding concentration (Fig. 4) . The upper boundary 
for the seeding ooncentration being set hy a vedue above vhich a speckle 
pattern is formed due to the interference of light scattered by the individual 
particles. If C is the particle ocxxrentration and AZ the width of the laser 
sheet, we obtain for the maximun ooncentration 

» 6p (6) 

where dp is the peirticle diameter. Qonsidering typicad vedues for dp = 0.01nm 
and Az = 0.5iim we obtain Cp<<5.10^^ m”^. 

The lOdier end for the seeding ooncentration can be determined by 
considering that, in order to have a valid measurement, a minimun nutiber of 
particle image pairs must be present in the area scanned by interrogation 
beam. The case of having a single particle image pair in the interrogation 
area is an ideal one, because it yields fringes with optimun SNR (Fig. 2). 
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However, this situation can only be achieved by lightly seeding the flow, thus 
giving rise to the so-called drop-out problem. An interesting case occurs 
vhen only two particle pairs are present in the interrogation area. The 
corresponding diffraction pattern includes multiple equally intense frixige 
patterns due to cross interference of non-corresponding image pairs. In this 
case, the local di^laceni^t cannot be resolved (Fig. 5). Finally, as the 
noiiber of particle image pairs in the interrogation area increases, the cross 
interferaice fringes become weaher in oonparison with the "main fringe 
pattern" vhich reflects the local displacanait (Fig. 6). The 
cross-interference fringes are sometimes designated as a "background speckle 
noise." CXir experience shows that, for reasonable fringe quality, at least 
four particle image pairs need to be present in the interrogation area. Ihis 
nuniber can be scrnevAiat relaxed if more than tvgo exposures were used for the 
photography. Considering that these particles exist in Idle volume defined dy 
the width, AZ, of the laser sheet, scanned by the beam vtose diameter is D, 
we obtain for the average concentration: 


Cp 


16 
AZ TT 




3.2 Exposure Parameters 


(7) 


The exposure parameters are chosen in accordance with the maxiimm 
expected velocity in the field and the required spatiad resolution. The 
spatiad resolution, vhich is equal to the cross sectional area, of the 
interrogating laser beam, is dictated by the scales involved in the fluid 
motion and should be less thaui the anallest scale. 

The time between exposures, T, is determined by the maximum permissible 
displaceraait of a peurticle such that a correlation is obtained vhen lcx:adly 
anadyzing the film negative with the probe laser beam. A necessary condition 
to obtain Young’s fringes requires that the distaunce between adjacent peurticle 
images be less than a fraction of the analyzing beam diameter. In practice 
the maximum permissible displacement that can be detected corresponds to the 
case vhen the fringe spacing, dp, is laurger than the diffraction limited spot 
size, d , of the interrogating optics. In analytical terms: 

/- > a - -A. 

where M is the magnification factor, D and X are respectively the 
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interrogation beem dieroeter and laser light venrelength. Vte obtain for the 
time betMsen exposures 


T = 


0.5 D 


Mv 


max 



( 9 ) 


is the maxlmun esqjected velocity in the field, and F the frequency 
oenrresponding to the inverse of the period T. For practical purposes v« use 
the ocxistant 0.5 instead of 0.8 as it vould be given by equation 9. 

Ihe exposure time, t, is a free parameter vdien a Gfj laser source is used 
together vdth a Gutter. In previous applications (3-7) the exposure time has 
been kept as anaiLl as possible for minimua peurticle image size. In a later 
study, Lourenco (11) demcxistrated that the duration of tl^e exposure time can 
be of the same order of magnitude as the time between exposures. 

For very short exposures the recorded particle images are identical to 
the diffraction limited particle images as the particles appear to be 
stationary during the ei^sure. For longer exposures the recorded images 
beoone streaks vAnose length is directly proportioned, to the exposure time. 
Ooherent illunination of these images generates Young's fringes superposed on 
a diffraction pattern which depends on the stree»k length. In the limiting 
case of a very short exposure, the diffraction pattern is syimetric and has a 
circuleu: shape (Fig. 2) . When the exposure time is increased, aimll streaks 
are generated, and the diffracted light in the spectrum is concentrated in a 
band vhose width is inversely proportional to the streak length (Fig. 7) . 

Because of the finite width of the diffraction pattern, the number of 

fringes, n^., in the pattern beoones a function of the ratio between the streak 

separation to the streak length, 1 . In analytical terms: 

s 

Op » + 1 (10) 

For an accuurate measurement of the fringe spacing, it is necessary to have a 
reasonable nunber of bright fringes. Usually five bright fringes give good 
results. 

Another parameter determining the degree of accuiracy with vhich the 
fringe spacing cam be determined is the contrast cind sharpness of the Young's 
fringes. If the number of exposures, N, is larger than two and the particles 
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the sharpness of the 


are displaced by the same anount after eadh exposure, 
interference fringes can be greatly improved. 

Par a multiple ejqxssure, the function describing the light intensity of 
the point sources beocmes 


N-1 

D(x,y) A *(^»Y + nMiy) (11) 

n=0 

resulting in an amplitude distributim of the transmittance of the negative in 
the Pxiri«: plane of a leis, as follows: 

^ ^ { N-1 K 

|D(u,v)|"‘ ^N+2 T] noosC(N-n)2 irMdy/XT) (12) 

' n»l ^ ' 

For a number of exposures, N, equal to two a relation similar to equaticm 4 is 
obtained expected. 

Figure 8 displays the fringe vd.dth for multiple exposures, normalized 
with the oorrespbnding width for a double exposure, versus the nuiber of 
exposures. It is clearly seen that an pptiimxn fringe width can be obtained 
with a limited nuiber of exposures. 

Ihe sharpening of the fringes using a multiple exposure recxsrding 
technique is of great eKivantage because it edlows a very easy visu 2 d and 
qualitative "measurement" of the direction and spacing of the young's fringes. 
It also increases their signed to noise ratio. Note that this improvanent in 
the signal to noise ratio is a result of the artificial increase in particle 
concentration due to the multiple exposure. 

3.3 Choice of the Film and laser Power Requirements 

The technique relies on the ability to detect and record on photographic 
media the seeding particles images. The particle is a function of the 
scattering power of the particles within the fluid, the amount of light in the 
illuninating sheet, camera lens and film sensitivity at the wavelength of the 
illvininating laser light. Although the pazt.icle detection increases 
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proportionally with increasing power of the illuninating las^, it is of great 
in(x>rtance to keep the laser power requirement to its minimum. An irportant 
reason for this limitaticHi is economy as the main cotponent in the cost of the 
apparatus is the laser source and its price increases very rapidly with the 
power delivered. 

In this section, we suggest films to be used and the laser pcwer required 
to expose these films. 

lb conpensate for the limited illuminating laser power, the films used 
to record the particle images, are required to have good sensitivity, but 
without sacrificing film resolution. When good precision is necessary, films 
with a resolution of about 300 line-pairs/nni and sensitivity from 25 to 125 
ASA (Afga Ortho 25, Kodak Technical Pan 2415) are suitable for the use either 
with pulsed Rtby lasers ( X~700rin) or pulsed, frequency doubled, NdYag lasers 
and CW Argon-Ion lasers ( A~500nn) . When a NdYag or a CW Argcxv-Ion laser are 
used more sensitive films may be used. A peurticular film, Kod£dc Royal-X Pan, 
1250 ASA, has excellent sensitivity in 400-€00nii range, and retains fairly 
good resolution (100 lin&-pairs/nm) . 

Another parameter to be accounted for is the film grain. When 
illuninating the film negative to produce Young's fringes, the film’s 
unexposed grains can introduce amplitude and/or fhase chcinges into the 
wavefrmt of the analyzing laser beam and thus create additional noise. This 
noise has genereilly a frequency content which is in the same range as the 
Young's frir^es. Therefore its elimination is difficult either fexin opticed 
or digital filtering. Film grain Ccuinot be totally eliminated and it will 
always be present in different degrees, especially in applications where fast 
film is used, to cope with low power density of the illuminating sheet. A 
method of soluticxi as been proposed (6) by which this source of noise is 
considerably reduced. It consists of producing a positive copy by contact 
printing on a very high resolution fine grain film, for example holographic 
plates. The positive copy is analyzed in the same manner using the laser 
probe beam. Figures 9a) and 9b) display resp>ectively the interference fringe 
piattems generated for the film negative and corresponding positive. A 
remarkable decrease in noise level is achieved. Another adveuitage of using 
this psrocedure is that the intensity of central zero order spot, due to the 
laser pardbe beam, is also decreased or eliminated. 

Once the film is selected, the amount of laser power required for the 
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exposure can be caiputed. For a successful exposure of the photographic 
emulsion the energy scattered by a tracer particle and vieved through the 
camera lens has to be larger than the film sensitivity at the vavelength of 
the llluninating laiser. It is custanary to write this condition in the 
following analytical fom: 


E 


■/’ 

“Tit 


Tdt>CE_ 


(14) 


Where E is the mean exposure level for a single particle (j/m?)^ T ig 

average intensity of the light scattered by a particle (W/m ) , E 0 is the film 
fog level and At is the ej^sure time. C is a constant between 1 and 10. The 
fog level is defined as the exposure level bela/ v^ch the transmissivity of 
the film is independent of the incident intensity (Fig. 10). Ihe average 

^ - T- ri'f 1 TrrVi+’ Vwr a anrw^lA ^ A 

terms of ( 10 ) 


T 



/• 


dft 


(15) 


2 

v^re Iq is the intensity of the illuninating ^eet (W/m ) , k is the 
illunlnating laser light wavenunriOer, a - o • a * is the Mie parameter, the 
vievmngle of the cam^a lens and d^ is the dimension of the diffraction 
limited particle image ocmposed of tvo terms ( 10 ), neglecting aberrations and 
grain noise, as follows 


=V'*p * 

Ihe first cxie being the particle diameter and the second one the edge spread 
caused by the limited respoise of the recording qptics 


dg = 2.44 X(14M)F# (17) 

with F# the camera objective F nvmber. 

When a pulsed laser is used the laser power required is sinply determined 
frcm equation 14, considering that the particle is stationary during the 
exposure 
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(18) 


Io> 


CEq( irK^d^) 

Yo2 dfi” 


Considering that for a pulsed laser tihe pulse duration (» exposure time) is a 
fixed peuraroeter, it is custanary to express equation 18 in terms of energy per 
pulse, 


e > 


o 




(19) 


Uie recxmnendecl value for the ocaistant C is between 3 auid 5 to ooipensate for 
film reciprocity effects caused the very short exposure. 

In applications involving the use of CW laser the exposure time is 
variable and considerably longer. In this case integration of equation 14 is 
saneMhat more eledxsrate. For the sake of sinplicity let us assune that the 
fluid motion, and hence the seeding particles displacement is unidimensiaial . 

r^, r and u be respectively the initiad position, current position and 
velocity of a tracer particle, in the film plane. Ihe current particle 
position can be expressed as 


r = r + u t 
o 


( 20 ) 


The exposure can be determined integrating equation 14 and oonsidering two 

d. d. 

cases, t < and t > ~ 
u — u 



T 


u 


T At 

I ( 


T 

T 

T 

u 


u 

d. 

1 

u 

(r -r+ 2d. + u 
o 1 


At) 


r <r<r +u. At 
0—0 

r +uAt<r<r +d. 

O — O 1 

r 4d. < r < r +d. + u. At 
o 1 — o 1 

( 21 ) 

O — O 1 

r +d < r < r +d.+ u. At 
o — O 1 

r +d4-u.At < r < r +2d. +u.At 
o — o 1 
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d. 

*Rie tnaxiimni eaqpoeure is obtained foe expoavre tiroes larger or equal to 

Thus, for the vAiole field, the optiroun exposure, is equal to d^/V ^^ vbich 
results in the maximum exposure, %d.thout a significant reduction on the 
dynamic range of the technique, as shcMn in the following section. 

Hie Of l 2 iser power required is given by. 


“o'" 


K^a?)v 


0-4 


/” 


max 




da 


( 22 ) 


Equation 19 diows that OH laser power required for a measurement increases 
with the flow velocity. Hiis limitation does not apply to the case of a 
pulsed laser vhere the pulse duration is independent of the flow velocity. 
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4 . VB[£cnY msnac baixz 


In the following we will analyze the technique's abilities of resolving 
large velocity gradients in flow fields, i.e, its dynamic range. The dynamic 
range is m the largest velocity difference that can be detected in 

the flow field. 

The low end of the dynamic range is determined by considering that for a 
measurement the spacing between the successive particle diffraction limited 
images or streaks is well resolved, i.e. do not collapse on each other. In 
anadytical form 


1 » d. + V.t < V.T 
s 1 


(23) 


Oonsidering that the time between exposures, T, is a function of the maximvam 


velocity substituting t - o for a pulsed laser and, t = 


di 
V 


for a CW 


max 


laser, yields the following expressions for the minimun velocity, 

2 M V d. 

V . =* ' pulsed laser 


nun 


(24) 


min 


2 M V d. 
max 1 

D - 2 M d^ 


OV laser 


(25) 


The velocity difference AV = — ^ — 

''min 

obtain. 


V _ V . 

max - velocity dynamic range and we 


AV = 


D 


2 M d^ 


- 1 


D - 2 M d. 


pulsed laser 


CW laser 


(26) 


(27) 


The dynamic range increases with shorter pulse durations and is maximun in 
applications involving the use of a pulsed laser. Considering typiczd values 
for d^ = .03nro, D = 0.5ram and M = 1, we cbtain for the dynamic range, 

AV 2. 7.5 pulsed laser (28) 

AV ~ 6.5 CW laser (29) 
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5. NJoamcf OF THE mmiciuE 


Hie overzdl accuracy of the technique depends on the accuracies vhich can 
be achieved in ^le photographic procedure and processing techniques. In this 
report we only discuss the accuracy of the photographic recording. The 
sources of error are the lens aberrations and the limited film resolution, 
causing the position of the particles to be recorded with an inherent error, 
and the spurious oc^trlbutions on the in-plane displacement recording the 
out-of-plane motion. 

In section 3.1 it vas pointed out that out-of-plane motion was a severe 
limitatim on the use of spiecicle in Fluid Dynamics applications. Ihe reason 
for this limitation was that slight out-of-plane motion by the scatterers, 
between the multiple exposure, results in non-identical shifted patterns 
poorly correlated. To cope with tiiis problem cne cerates in the Particle 
Imaging mode, vhere the particles are directly imaged and recorded. However, 
as shown in the following, out-of-plane motion may also ocxitribute to 
considerable errors in the velocity measurement. 

Let us consider the imaging system of figure 11 emd the particle in 
positicxi within the laser sheet. Between exposures the particle moves to a 
position M^, dxie to three-dimaisioned. fluid motion. Ihe oc mp onents of the 
displacement vector are dx, dy, dz, with dx and dy, the in-plane and dz the 
out-of-plane ocn^nents of the displacement. In the image plane the 
photographic plate records the pxjsition of the particle at Pj^ and The 
coordinates of these points are given by: 


^»!(x+dx) (1+ dz/d^) 
-M(y+dy) (1+ dz/d^) 


^ 6 ^ 




vihere second order terms have been neglected for siitplicity. 

Ihe displacemoit measured by means of Young's fringes is given by 


<l * ) 


ay^ = Mdy (1 + ^ ) 


eind the measured displacanent ocmponents referred to the object plane are 
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*-0 


dx + 


xdz 


■% - dy + ^ 


The oontribut.ion of the out-of-plane displacement to the measured 
displacement is given by two parasite terras (xdz/dj^) and (ydz/d^^) . The error 
produced by the out-of-plane rooticxi while negligible in the neighborhood of 
the optical axis, increases linearly, and may beccrae inportant, with the 
distance fron the c^iced cucis. The influence of the out-of-plane motion 
beoctnes particulcurly in^rtant vhen imaging the flow with short focal and wide 
angle cbjectives. 

A relation can easily be found for the theoretical error in a measurement 
due to out-of-plane notion: 




(33) 


vhere ip and ij/ are respectively semi- field angle oc mp onents alcxw the x 

A y 

and y axis respectively. Hence the error depends only on the ratio of the 
out-of-plane to in-plane ocnponent of displacement and the tangent of the 
semi- field angle. Figure 12 shows the theoreticail relative error for various 
ratios of dz/dx euid tam ipx. 
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6. VALUKTIQN OP TIE TBCaUQUE 


6.1 Exp^ijnental Oonfiguration 

In order to evaluate the present capabilities of the technique a bench- top 
experiinent %<as perfoimed. In view of the application of the technique to 
vortical flows, two flows that include inportant vorticeil motions were 
generated and recorded. Both flows were created by towing models in the 
reduced scale Fluid Mechanics Research Laboratory towing tank facility (Fig. 
13). H»e first cxie consisted of the Kazman vortex street generated by a 
cylinder with 10 nin in diameter, and the second one by the flew over a 60nm 
(dK>zd airfoil (NAOV 0012) , at a 30 degree incidence. Both models were towed 
with a velocity of 23.5 mn/sec. The fluid used in these experiments was vater 
seeded with small, .004 mn in diameter, particles (TSI model 10087). The 
oorrespcxiding Reynolds nunbers were 230 for the cylinder and 1400 for the 
airfoil. These flows are excellent test cases because they include large 
scale vortical motions and extreme velocity gradients. These extreme 
ngpradimits serve as a test to the technique's capeibilities of providing 
information over a large velocity range. 

A view of the laser sheet arramgement is also shown in figure 13. The 
laser beam fron a 5 Vfeitt Argon-Ion laser (Spectrar-Physics series 2000) is 
steered and focused to a diameter of .3nm using an inverse telescope lais 
eurrangement. A cylindriced. lens, with a focal length of -6.35 nm, is used to 
diverge the focused beam in one dimension, creating a light sheet. The laser 
sheet is 70 ran wide and illuninates the mid-span section of the models. 

For the multiple exposure, the CW laser beam is modulated using a Bragg 
cell. This shutter has been assembled using stamdard equipment. A schematic 
is shown on figure 14. The advantages in using this shutter arrangenoent are 
the possibility of an ind^ndent cihoice of the number of exposures, the time 
between exposures and the exposure time. 

Due to conditioning cptics and Bragg cell losses, the available laser 
pamec was reduced to 70% of the total power. In this experiment the laser 
delivered 8 V^atts in the multiline mode. The power density, I^, of the laser 
sheet was 


'o “ ' .27(wW) 


(34) 
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with d|^ bZ the laser sheet czx^ss section, vAiere » .3 inn is the width of the 
focused leiser beam, and AZ » 70 inn the laser sheet span. 

The recording optics ocnsisted of a 4x5" format cannera with 
interchangeable objective lens. In this experiment we used a 135 mn Imis. 

The lens aperture vias set at F# 5.6 and the magnification factor was 0.47. 

The cemera was positioned in the ledxxratory referaice frame. Therefore, the 
flow field was recorded in the reference freme moving with the mean flow. The 
film used was a N3EAK techniced. pan 2415 with a sensitivity of 125 ASA. This 
film has a resolving power of 320 line-pairs/mn and extremely fine grain. The 
large resolving power is required for increased recording accuracy. This 
corresponds to a inaccuracy of 0.002-^.003 inn in the recording of the position 
of the tracers which corresponds in the present case to an error of 1% or less 
in the high end of the velocity range. 

6.2 Results and Discussion 


6.2.1 Flow behind a circular cylinder 


The technique was first tested in the measuranent of the flow bdiind a 
cylinder. The frequency of exposures was optimized according to equation 9 
and equal to 90 Hz. The exposure time was 1 msec, which corresponded roughly 

Figure 15 is a triple exposed fhotograph of the flow. Characteristic 
fringe patterns, resulting from coherent illunination of selected locations of 
the photograph are also displayed. 


The velocity data was acquired in a regular mesh by digital processing of 
the Young's fringes, produced by point by point scanning of the photograph. 
The scanning step size and Idie dimension of the analyzing beam diameter were 
.5 nm, which corresponds to a spatial resolution of about 1 ran in the object 
plane. The fringe patterns wmre processed using the interactive 
one-dimensional averaging software descrfoed in references 4,6 at the von 
Karman Institute Digital image Processing Ehcility. 

Figure 16 shows the mapped two-dimensional velocity vector field 
superposed on the original photogreiph. As it can be observed, this data 
represents with great fidelity the flow field. However, in sane regions 
velocity data drop-out occurs, due to the following reasons. Firstly, at the 
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lcx:!ation8 near the instantaneous centers of rotation of the flow, the velocity 
of the particles is very close to zero and beyond the low end of the dynamic 
raur^e of the technique. This could be avoided if the convective flow velocity 
vas measured, for exanple, 1:^ having the camera traveling with the model. 
Seocxidly, .some data drop-out occurs at random locations due to various reasons 
such as scratches in the film, fringes with bad SNR or lack of seeding. These 
regixxis of signed drpp-out were "filledUn" using an interpolatim algoritlm. 
Because the two-dimaisional velocity field is measured on a squaure mesh, the 
spatial derivatives of the velocity, and consequently the vorticity, can be 
estimated by a sinple algorithm. Oonsidering that the grid location of the 
meatsureroent points are labeled with indices i and j, the vorticity con^nent, 
at location i, j is 


“ i'j 


_ ^itl,j ~ ^i~l, j ^i>i-H ~ ^i,i-l 
2 Ax 2Ay 


(35) 


For the boundaries of the velocity field, an excentered scheme is used to 
evaluate the spatiad derivatives. This results in a smadl inatccuracy on the 
evaluation of the vorticity at the boundary. Figure 17 shows the vorticity 
contours for the flow behind the cylinder, liiproved visualization of the 
velocity and vorticity fields is obtained by color encoding the different 
levels of the vorticity ocnponent ft as shown in figure 18. 


6.2.2 Flow over the iJACA 0012 Airfoil 


For the flow over the airfoil the frequency of exposures vas varied over a 
wide range from 60 Hz to 90 Hz. Preliminary analysis of the film negatives 
showed that best results were obtained with a frequency of 90 Hz. The 
exposure time was 1 msec and kept constant throughout the experiment. 
Ehotografhs obtained with different nunber of exposures eure dXMn in figure 
19. The figure demonstrates that a significant increase in quality of the 
'•flow visualization" is achieved with a larger nunber of exposures. Another 
equally important effect of the nufoer of exposures is the iirproveraent of 
fringe sharpness and SNR. 

A ocnp»lete analysis of the triple exposed photograph (Figure 19 b) ) vas 
also carried out and results are presented on figure 20. 
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7. OGKUBIClfi 


Using tile apparatus described in the previous sections a series of 
Particle Displacement Lnages were successfully taken usii^ different exposure 
paremeters. local coherent illunination by the probe leiser beam yielded 
Young's fringes of good quality at almost every location of the flow field. 
Using the von Karman Institute image Processing facility, these fringes were 
analyzed and the velocity and vorticity fields were derived. Ihe present 
stuiy has lead to the following oonclusions: 

(i) Particle Image Displacement Velocimetry is inieed suitable for the stuiy 
of flows with vorticed motions. 

(ii) In the present study, the flow velocities were limited to a anzdl value 
(50nn\/sec) because of laser power limitations. However this limitation 
does not apply if a pulsed laser delivering high energy pulses is used. 
CXx:e a pulsed laser becomes available to us, tests in the IMRL wind 
tunnel will be conducted to prove the feasibility of the technique in 
higher velocity (up to ItHSm/sec) air flows. 

(iii) A larger towing tank facility has just become available to the BMRL in 
which tests oan be carried out in a wide range of towing velocities. 

This facility is provided with camera mounts that travel with the towed 
model. Therefore it is possible to capture the flow in two reference 
frames: the one traveling with the model and that of the laboratory. 

(iv) The next step in the development of the technique will consist on the 
development of schemes to ocnpensate for the error introduced by the 
out-of-pleme motion and device methods for its correction. 

Simultaneously a method for the measurement of the magnitude of the 
out-of-plane motion will be developed. 

(v) The PMRL cap^ilities in the area of digital processing of the Young's 
fringes is paresently being developed. In a later stage Image 
Ootrelation techniques, will also be investigated, and a comparison of 
the p>erformance of the two methods will be ceoxied out. 
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Figure 1. Schematic of 
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Figure 3. Schematic of the image correlation function 
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Figure 4. Mode of operation versus seeding concentration 
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Figure 6. Fringe pattern generated by multiple image pairs. 
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Figure 8. The variation of fringe width with number of exposures. 



Figure 


Young's Fringes 



a) film negative b) film positive 


9. Young's fringe patterns obtained with film negative 
and positive. 
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Figure 10. The density- exposure curve for a film 
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Image Plane 


Figure 11. A schematic of the imaging system for the 
estimation of the out-of -plane motion. 
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Figure 12. Theoretical re 
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Figure 13. Photograph of the experimental facility. 
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Figure 14. Schematic of the laser shutter 
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Figure 15. Triple exposed photograph of the vortex street. 
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Figure 16. Velocity vector map of the vortex street 
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Figure 18. Velocity and vorticity fields of the Karman vortex street. 
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Figure 19c. Multiple exposure photographs of flow past an airfoil 
at 30 degrees angle of incidence; a) double exposure, 
b) triple exposure, c) Quadruple exposure. 
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Figure 20. Velocity and vorticity fields of flow over a NACA 0012 airfoil 




